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We demonstrate experimentally a precise realization of Coulomb Blockade Thermometry (CBT) 
working at temperatures up to 60 K. Advances in nano fabrication methods using electron beam 
lithography allow us to fabricate a uniform arrays of sufficiently small tunnel junctions to guarantee 
an overall temperature reading uncertainty of about 1% 


INTRODUCTION 


and for the limit Eq ^ ksT laii] with 


Primary thermometers are required for the ongoing 
effort to redefine the Kelvin scale based directly on a 
fundamental physical constant, namely a fixed value 
of the Boltzmann constant /cb [I]. This task requires 
primary thermometer realizations using different ap¬ 
proaches, where consistency between different methods 
is compulsory. Noise thermometers, like the Magnetic 
Field Fluctuation Thermometer (MFFT)[2l [3], Current 
Sensing Noise Thermometry (CSNT) [4] and the CBT [5] 
are actually candidates that are under consideration for 
the low temperature range. Experimentally, the method 
of CBT is a well established tool for precise determina¬ 
tion of temperatures below 1 K [5j [6] with good precision 
[7]. At higher temperatures though, fabrication inhomo¬ 
geneities [8] are setting strict limits to the feasibility of 
CBT, so that no precise sensors operating above about 
10 K were reported so far. We will present in this work 
a CBT operating at temperatures up to 60 K and dis¬ 
cuss the uncertainty component arising from fabrication 
imperfections. 


THEORETICAL BACKGROUND 


Ui /2 = b.439NkBT{l + 0.3921^)/e, (2) 

Ctt 

In addition, a secondary temperature measurement 
can be obtained by recording the depth of the zero bias 
conductance dip 
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with the defined parameter un = EclksT. Equation (§ 
has to be extended by higher order contributions when 
leaving the strict limit of Ec ksT. The measured 
conductance follows well the expression [6] 
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We use in this work a numerical [5] way to precisely 
calculate the conductance curves that includes no ap¬ 
proximation at any ratio of Ec to ksT in order to extend 
the description into the intermediate Coulomb blockade 
regime {Eq — ksT) [6]. This allows us to calculate 
conductance curves also including the overheating effects 

m- 


Coulomb Blockade Thermometry employs the temper¬ 
ature dependent electrical conductance G of tunnel junc¬ 
tion arrays. Classically, CBT works in a weak Coulomb 
blockade regime Eq ^ ksT, where the charging energy 
of the system with N junctions in series is defined as 
Eq = [(A^ — 1)/A]e^/G, where G is the capacitance of 
one island [5] . The latter depends on the physical size of 
the tunnel junction in combination with the thickness of 
the insulator, in addition to the stray capacitance of the 
island. Note that temperature sensors for an increased 
temperature range need to have an increased charging 
energy and therefore smaller junctions. Temperature is 
derived in the strict limit of Eq ksT from the mea¬ 
surement of the full width of the conductance dip at half 
minimum using 

( 1 ) 


FABRICATION OF THE SUSPENDED 
GERMANIUM MASK 

We fabricate the devices using electron beam lithogra¬ 
phy combined with shadow evaporation technique pT[[T2] 
and a special tri-layer resist scheme [Hill] that employs 
in our case a 22 nm thin germanium (Ge) film as the 
material for the suspended mask. The wafer preparation 
is done as follows: at first, we spin coat (5500 rpm for 
60 s) a 4 inch silicon wafer covered with 300 nm of sil¬ 
icon oxide with a 400 nm thick copolymer film (11 % 
MMA/PMMA in ethyl lactate) and bake it for 45 min 
at 175 degree Celsius. This extremely long and hotter 
than usual baking time is required to avoid out-gassing 
of the copolymer underneath the gas-tight Ge layer that 
is deposited next. We use an electron gun evaporator 
equipped with a powerful cryo-pump to deposit the 22 


= 3.A39NkBT/e. 


2 


nm thin Ge film with a growth rate of 0.1 A/s in order to 
achieve a low stress Ge film. This corresponds to a total 
evaporation time of about 22 minutes for the 22 nm film. 
We have to control the deposition rate precisely since it 
has a strong effect on the strain of the Ge film, which 
can vary from positive to negative. The optimal rate is 
stable in our evaporator, but in general depends on the 
vacuum conditions. We observe that more perfect vac¬ 
uum conditions require a lower deposition rate. Finally, 
we again spin coat the wafer with a 50 to 60 nm thick 
standard PMMA electron beam resist (molecular weight 
2200k) using a 2.25 % solution in anisole at 2500 rpm and 
bake at 150 degree Celsius for 90 seconds. This proce¬ 
dure yields a very reproducible quality of the resist stack 
and allows us to achieve a constant resolution and device 
quality over the full 4 inch wafer surface. We perform the 
electron beam lithography with an e-beam writer having 
100 keV electron energy, doses of 1000 /rC/cm^, no prox¬ 
imity correction and beam currents of typically 1 nA. 

The development for the mask is separated into three 
phases [131 |T4]: firstly, we develop the exposed top layer of 
PMMA using standard MIBKiIPA 1:3 solution, followed 
by pure IPA bath and finally rinsing in deionized water, 
30 seconds for each process step. Secondly, we transfer 
the opening in the PMMA to the germanium film using 
reactive ion etching (RIE) with CF 4 plasma at a pressure 
of 40 mTorr and a rf-power of 40 W for 140 s. The CF 4 
plasma selectively etches only the germanium film and 
does not affect the resist noticeably. Thirdly, we use 
pure oxygen plasma at a pressure of 30 mTorr and a 
rf-power of 40 W for 15 min that an-isotropically etches 
selectively the copolymer film and removes at the same 
time the PMMA film covering the remaining germanium 
mask. Increasing the pressure of the oxygen plasma to 
225 mTorr for additional 20 min etches more isotropically 
and creates the necessary undercut below the germanium 
mask, enabling shadow evaporation at an angle of about 
± 25°. 

The next step in the fabrication is the shadow evapora¬ 
tion performed in an electron gun evaporator: we deposit 
first under a tilt angle of + 20 ° a 20 nm thin aluminium 
film with a high deposition rate of about 12 A/s. The 
high deposition rate produces a film with a reduced grain 
size. Then, we oxidise the aluminium film with 5 mbar of 
pure oxygen in the deposition chamber for 300 s to create 
the AlOx tunnel barriers. Finally, we deposit the second 
aluminium film with a tilt angle of -20°. The lift-off is 
done by immersing the wafer in acetone at ambient tem¬ 
perature. Figure depicts a small section of one sensor 
that consists in total of an array of 77 junctions in series 
to increase the signal strength and to suppress the influ¬ 
ence of the electromagnetic environment [H [15] and in 30 
parallel rows to reduce the final resistance of the sensor to 
a value of about 100 kfl in order to facilitate the electrical 
measurement. Already the design of the mask takes the 
main requirement of the sensor, the small junction size 



FIG. 1. A scanning electron micrograph of a sample identical 
to the measured one. Depicted are three rows of 6 junctions 
each fabricated at two angle shadow evaporation (see text). 
We observe a remarkable reproducibility over the full 77x30 
junction array, note that the 2nd shadow, that falls in the 
image always on the left, has a slightly reduced size due to 
partial blocking of mask from the first angle evaporation.) 

in combination with a small scatter in the junction areas 
into account: the crossing geometry, where the second 
shadow overlaps the first one totally defines the junction 
area leaving only the line width fluctuations as an error 
source. Even smaller junction areas are possible to fab¬ 
ricate if the second shadow would only partly cover the 
first one, but in that case, both the length of the line and 
a misalignment of the mask with respect to the tilting 
angle would contribute to the error in defining the junc¬ 
tion size and consequently the important tunnel junction 
resistance. As depicted in Fig. the final junction size 
of about 45 nm x 32 nm is constant over all sensor junc¬ 
tions, the material deposited for the first layer blocks 
partly the mask and yields the reduced line-width of the 
second shadow. 


NONLINEAR BACKGROUND CORRECTION 

FOR THE OPERATION OF CBT SENSORS AT 
ELEVATED TEMPERATURES 

The measured total sensor resistance is 78 kO (see Fig. 
1 ^ corresponding to a resistance per junction of 200 kO. 
A practical realization of a CBT sensor shall show a linear 
voltage dependence up to applied bias voltages that are 
ideally few times the half width of the Coulomb peak. At 
low enough temperatures, the assumption of a bias volt¬ 
age independent conductance (apart from the Coulomb 
blockade peak) of the aluminium oxide tunnel barriers 
holds well enough: the Coulomb peak is well separated 
from the non-linear background even at a temperature 
of 8 K (see Fig. [^. 

Real devices show due to the finite tunnel barrier 
height (about 2 eV [T 6 | in the case of a tunnel barrier 
made from AlOx) noticeable deviations at higher volt- 
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FIG. 2. (left) Measured voltage dependent differential resis¬ 
tance curves at two different temperatures for a CBT sensor. 
The curve measured at 8 K was clipped to emphasize the non¬ 
linear background, the zero bias value reaches 95 kQ. Both 
uncorrected and corrected curves are depicted, (right) Zoom 
into the area affected by non-linear background, the blue line 
is the function describing the background and the black hor¬ 
izontal line is a guide to the eye at the position of the tunnel 
junction resistance for small bias values. 


ages: such a finite barrier height leads according to the 
phenomenological Simmons model m to a parabolic de¬ 
viation in the form G{V) = Go(l + {V/Vo)‘^) with a typ¬ 
ical value of Vb = 0- 25r[ini[IH]. Aiming at a temper¬ 
ature reading of 60 K requires a bias voltage range of 
i^BiAS ~ 2 Vi/2 ~ 50 mV per junction as increasing the 
temperature also increases linearly V 1 / 2 , where the sim¬ 
ple model predicts a noticeable deviation on the order 
of 1%, whereas the peak-hight itself is reduced to about 
3% in the sensors discussed here. Experimental findings 
differ somewhat from the Simmons model, even a more 
sophisticated model [19] employing first principle calcu¬ 
lations to describe the density of states in the tunnel 
barrier catch well the background at high voltage, but 
falls short in describing the detailed behaviour close to 
the region of interest here. As a solution, we determine 
the background using the depicted curve measured at 8 
K (see Fig. by a polynomial function in the form of 
R = Ro{I (see Fig. (right)) and correct the mea¬ 
sured curves with this function at all temperatures. An 
alternative way to limit the fitting to the range where the 
background is flat (+/- 30 mV/junction) reveals within 
the confidence interval of the fitting routine no deviation 
between the two methods. Figure (right) shows, that 
the course of the background is very similar for differ¬ 
ent temperatures and the applied correction brings both 
curves to a constant value within the noise level of the 
measurement. In conclusion, we show in this section that 
we can neglect the uncertainty component arising from 
the background at low enough temperatures 40 K), as 
the background does not deviate noticably in this range. 


Furthermore, the presented background correction is ef¬ 
fective at temperatures between 40 K and 60 K. 

EXPERIMENTAL CHARACTERIZATION OF 

THE MULTI TUNNEL JUNCTION SENSORS 

We use a standard setup [ 7 ] including lock-in amplifiers 
to measure the bias dependent differential conductance 
curves of the devices. The sensors are wire bonded to 
a sample holder that is mounted in the vacuum cham¬ 
ber of a dilution refrigerator immersed in a liquid helium 
dewar. The dilution refrigerator is operated without con¬ 
densed ^He/^He mixture in order to access the tempera¬ 
ture range up to 60 K: we apply heat to the sample stage 
with the standard electrical heater that overcomes the 
cooling towards 4 K due to the residual heat exchange 
gas in the vacuum chamber and the thermal conduction 
of the apparatus. 



FIG. 3. Measured (black dots) and numerical fits (lines, see 
text) of bias voltage dependent, normalized conduction values 
of a GBT sensor at four different temperature values between 
10 and 60 K. The voltage axis is scaled to depict the average 
voltage drop across one single junction in the sensor array. 
The top panel depicts the ratio between the measured (Gm) 
and the simulated (Gsim) conductance at 59.6 K. 

Fig. i depicts few measured conduction curves to¬ 
gether with the extracted temperature values from the fit. 
Despite the unusually high temperatures, the model still 
has to take overheating effects into account m as the 
volume of the metallic islands is quite small (about 5e- 
4 /i m^) and the applied voltage and consequently joule 
heating is increased: At a base temperature of 9 K, the 
electronic temperature in the islands rises at the applied 
voltage of Vi /2 by about 100 mK. Already the perfect 
agreement between the theoretical model and the mea¬ 
sured data, the magnitude of the signal even at exception¬ 
ally high temperature of 60 K and the flat background 
underline the outstanding quality of these devices. We 
find a charging energy as high as Ec/kB = 10.6 K, 
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meaning that the sensors operate even at this high tem¬ 
peratures mainly in the intermediate Coulomb blockade 
regime [6]. This is very beneficial for the signal magni¬ 
tude and for suppressing the remaining influence of the 
non-linear background. 

EXPERIMENTAL CHARACTERIZATION OF 
THE SENSOR UNIFORMITY 

The main topic of this chapter is the experimental 
characterization of the fabricated multi tunnel junction 
sensors with respect to the resistance uniformity be¬ 
tween the individual junctions. The latter is too difficult 
to probe directly because of the mere number of junc¬ 
tions and their small size. Instead, we follow here four 
approaches to determine this value experimentally and 
show that all yield a consistent result: (a) we fabricate 
together with the real device a test structure that enables 
to probe 11 individual test junctions, (b) measure three 
sensors on the very same chip in a range of temperatures 
to analyze the deviations, (c) compare the sensor to a 
reference temperature when immersed in liquid helium 
and finally (d) compare the temperature reading of the 
device to a single tunnel junction thermometer. 

Probing of individual junctions 

We performed a direct measurement of the tunnel junc¬ 
tion resistance at room temperature using the structure 
depicted in Fig. that is fabricated along with the real 
structures and allows monitoring of the fabrication accu¬ 
racy. The 12 leads connect to bonding pads that are wire 
bonded to the test set-up. 



FIG. 4. Scanning electron mikrograph of an array of 11 junc¬ 
tions connected to 12 leads allowing electrical probing of all 
junctions individually. The junction size and layout is identi¬ 
cal to the size of the structures in the sensors. 

The immediate environment close to the junctions on 
the chip during fabrication of the test junction array dif¬ 
fers somewhat from the real sensor (compare to Fig. Q 
due to the additional leads, but the inspection with a 


scanning electron microscope reveals no significant influ¬ 
ence on the junction area or quality of the mask. Con¬ 
sequently, we assume that we can use the results of the 
test structures to model the full sensor array: We find 
in the best case a mean resistance of 70 kU with a stan¬ 
dard deviation of 3.3 kU, corresponding to 4.7%. We 
use this result to simulate the deviation of a full sensor 
consisting of 77 junctions using the method described 
in [HI [20] and find a resulting temperature deviation of 
0.07%. As an upper limit, we found a resistance devia¬ 
tion of 11.4%, still resulting in a temperature deviation 
of 0.7%, well below 1%. These findings allow under the 
given assumptions to set an upper limit of the uncer¬ 
tainty component of temperature determination arising 
from fabrication non-uniformities to 1%. 


Inter—comparison of three sensors 

We describe in what follows the inter-comparison of 
three CBT sensors in the temperature range between 10 
K and 60 K. The result is at the first glance already very 
promising: the three sensors agree with each other on 
a level of about 0.2% over the investigated temperature 
range. 

The setup (Fig. [^left) is specially optimized in order 
to rule out other contributions to the measurement error: 
the very same single voltage relating the measurement to 
temperature is applied to all three sensors that are mea¬ 
sured at the same time. In addition, all three sensors are 
fabricated together with single junction devices on one 
single silicon wafer leaving only little room for tempera¬ 
ture gradients between the sensors. Finally, as the setup 
allows to measure the three sensors at the same time, 
any temperature fluctuations or drifts affect all sensors 
identically. This leaves only the error source arising from 
the sensors itself influencing the results, if one assumes 
in addition no significant deviations from linearity in the 
amplifier chain. 

We discuss two further questions at this point: (a) 
what is the expected absolute uncertainty that is consis¬ 
tent with the observed relative deviation? And (b) to 
what upper limit of junction array inhomogeneity does 
this value correspond to? In order to answer those ques¬ 
tions, we follow the approach described in m- we calcu¬ 
late numerically the conductance curves for arrays that 
have a random fluctuation of resistances in the individual 
junctions characterized by a Gaussian distribution with 
a set of different standard deviations. Now we can link 
the expected absolute deviation of temperature reading 
for a sensor to the relative deviation among many sen¬ 
sors with different individual random fluctuations of re¬ 
sistances. The result is that a 10% scatter in the junction 
resistance (corresponding to about 0.7% absolute devia¬ 
tion in temperature reading) would lead to a 0.2% devia¬ 
tion among different sensors and describe best the exper- 
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FIG. 5. (left) Experimental set-up measuring three CBT 
sensors simultaneously: one single voltage source is attached 
to three sensors. The voltage consists of two components, 
the resistance network adds a divided DC bias voltage set by 
Rg/Rdc and a sinusoidal AC component via Rg/Rac- The 
resulting differential conductance of the sensor is measured us¬ 
ing one current amplifier and one lock-in amplifier each. Only 
the sensors (blue shaded area) are located at low temperature 
and connected to room temperature via filtered Thermocoax 
lines, (right) Resulting relative deviation among three sen¬ 
sors (red, blue, green circles) with respect to the measured 
mean temperature value of three sensors in the temperature 
range between 10 and 60 K. Given exemplary error bars are 
the one sigma confidence interval of the fit to the measured 
conductance curves and agree well with the observed scatter 
of the temperature data. 

imental findings. These results are again consistent with 
an upper limit of the uncertainty of temperature deter¬ 
mination arising from fabrication non-uniformities to be 
smaller than 1%. 

Comparison of the sensor reading towards a 
reference temperature 
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FIG. 6. Voltage bias dependent conductance curve (dots) of 
a sensor immersed in liquid helium with a fit to the full model 
(red line). The green lines mark the intervals, where a simple 
polynomial fit of second order is used to determine the half 
width of the dip. Open circles mark the position of the half 
value of the conductance dip on the curve, the rhombus marks 
the conductance minimum (see text). 

perature using Eq. (© (4.385 K). Note, that the agree¬ 
ment between all values is within 0.5% and that there is 
an excellent match between measured data and the the¬ 
oretical curve. These results are again consistent with 
an upper limit of the uncertainty of temperature deter¬ 
mination arising from fabrication inhomogeneities to be 
smaller than 1% for the presented sensors. Note that 
the charging energy of the sensor presented in this chap¬ 
ter was lower (5.5 K) in order to adapt the conductance 
dip magnitude to the operation temperature close to 4.2 
K. Sensors having the higher charging energy discussed 
above would already suffer from the influence of random 
background charge distributions at this temperature [6]. 


We demonstrate in this section a comparison of the 
sensor reading to the reference temperature realized with 
a liquid helium bath in combination with measuring the 
vapor pressure of the cryo-liquid: the sensor is immersed 
in a pressurized liquid helium dewar at 1150 mbar, cor¬ 
responding to a temperature of 4.359 K according to the 
ITS-90 temperature scale. We neglect at this point tem¬ 
perature gradients within the liquid helium and the un¬ 
certainty arising from the pressure measurement (± 10 
mbar), as the latter corresponds to a temperature devia¬ 
tion of i 0.2 %. 

Temperature is extracted from the conductance curve 
measurement in two ways: fitting the numerically cal¬ 
culated conductance curve to the full data range (4.363 
K) and by determining temperature via the half width 
of the dip in the conductance curve using second order 
polynomial fits to the dip and around the two values neat 
half minimum as indicated in Fig. to determine tem- 


Experimental comparison to a single junction 
thermometer 

Finally, we discuss the comparison of the CBT sensor 
with a single junction device (SJT, [71 |T5l [21] ). The ob¬ 
vious advantage of the SJT over a multi junction device 
is that no scatter of junction resistance in the array in¬ 
fluences the result. Like this, the direct comparison of 
the CBT and SJT sensor shall reveal the influence of the 
non-uniformity. We choose an optimal temperature of 
about 9 K for the experiment, where on the one hand 
the overheating effects and on the other hand the influ¬ 
ence of the non-linear background are still negligible. 

An image of the SJT is depicted in Fig. left inset: 
one central junction is connected via four long arrays of 
tunnel junctions allowing a four probe measurement of 
the central junction and creating a suitable electromag¬ 
netic environment for precise thermometry m- We use 
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FIG. 7. Measured (black dots) and fitted (red line) bias volt¬ 
age dependent conductance curves for a CBT (left) and a 
single junction (right) at a temperature of about 9.3 K. The 
uncertainties of the temperature readings given in brackets 
are the one sigma confidence interval of the fit to the mea¬ 
sured conductance curves and neglect other components (see 
text). The voltage axis of the CBT measurement is scaled to 
that of a single junction. The inset of the SJT graph shows 
a scanning electron micrograph of a device that is similar to 
the measured one (see text for details). 

a standard setup [7] for the measurement of the differen¬ 
tial conductance; the main difference with respect to the 
CBT setup is a voltage amplifier that is needed to de¬ 
tect the applied bias voltage across the single junction. 
The main disadvantage, the small signal magnitude, of 
using only one junction is partly compensated by an in¬ 
creased time constant of the lock-in amplifiers, leading 
to extremely long measurement times for a single data 
point of about 12 hours in order to yield comparable 
signal to noise ratios for both sensor types. The latter 
causes some long term drift effects in the background of 
the SJT conductance curve (Fig. right), most likely 
due to amplifier gain deviations caused by room temper¬ 
ature changes. The agreement of the two temperature 
readings is perfect. One has to consider the main uncer¬ 
tainty components though: the voltage axis realization 
of the SJT relies on the accuracy of the amplifier gain 
that is limited to 1% whereas the main influence in the 
case of CBT comes from the sensor non-uniformity. This 
comparison can set again an upper limit of the accuracy 
for the described CBT sensors to 1%. 


the performance of the sensors due to the superior sig¬ 
nal to noise ratio, allowing for faster signal acquisition or 
higher precision. These advances have also consequences 
for sensors operating in the classical temperature range 
between 50 mK and 1 K: if one fabricates junctions with a 
larger area of about of 1x1 fivD? suitable for this tempera¬ 
ture range with the same process and absolute fabrication 
accuracy, the resulting variation of the total junction area 
will be negligible. Other effects, like variations in the ox¬ 
ide layer thickness [18] would become dominating in this 
case and would need to be optimized. 

We focus in this paper on two uncertainty components 
that are crucial for high-temperature measurements: the 
background correction, where we show that the uncer¬ 
tainty component can be made negligible due to an in¬ 
creased charging energy and the array uniformity where 
we demonstrate an upper limit of 1% in this work. The 
whole uncertainty budget was not assessed. Other impor¬ 
tant uncertainty components are the fitting uncertainty 
due to measurement noise and disturbances and finally 
the measurement of the dc voltage in the G(V) measure¬ 
ments. The latter can be suppressed roughly to 0.01% 

[23] . We can demonstrate in this work only an upper 
limit for the uncertainty, but our results would allow the 
conclusion that the presented sensors are even more ac¬ 
curate. A direct comparison to a realization of the ITS90 
temperature scale in combination with a traceable volt¬ 
age measurement scheme would allow a characterization 
with an one order of magnitude smaller uncertainty. 

CB Thermometry turns out to be a mature technol¬ 
ogy operating nowadays in a temperature range spanning 
over more than four decades in temperature from 60 K 
as shown in this work down to the low milli-kelvin range 

[24] . reaching electronic temperatures as low as 3.7 mK 
if the sensor is immersed directly in liquid helium of a 
dilution refrigerator [25|. 
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